Astronomy & Astrophysics manuscript no. herxlinf 170108r2 


© ESO 2008 


February 11,2008 





INTEGRAL observations of Her X-1 

D. Klochkov\ R. Staubert\ K. Postnov^, N. Shakura^, A. Santangelo^ S. Tsygankov^, A. Lutovinov^, 

I. Kreykenbohm' '^, and J. Wilms^ 



00 

o 
o 

(N 

X) 

(D 

oo 



O 



00 



O 

OO 

o 



Institut fiir Astronomie und Astrophysik, University of Tubingen, Sand 1, 72076 Tubingen, Germany 
Space Research Institute, Profsoyuznaya str. 84/32, 117997 Moscow, Russia 
Sternberg Astronomical Institute, Moscow University, 119992 Moscow, Russia 
INTEGRAL Science Data Centre, Chemin d'Ecogia, 16, 1290, Versoix, Switzerland 

Dr. Karl Remeis-Sternwarte, Astronomisches Institut, Universitat Erlangen-Numberg, Sternwartstr. 7, 96049 Bamberg, Germany 



Accepted 29/01/2008 



ABSTRACT 



Aims. We investigate the X-ray spectral and timing properties of the accreting X-ray pulsar Her X-I observed with the INTEGRAL 
satellite in July- August 2005. 

Methods. The data analyzed in this work cover a substantial part of one main-on state of the source. The short-time scale pulse 
period development is measured. X-ray pulse profiles for different energy ranges and time intervals are constructed. Pulse-averaged 
and pulse-phase resolved broad band X-ray spectra are studied. Spectral changes during X-ray dips are explored. 
Results. The X-ray pulse profiles are found to change significantly during the period of observations. For the first time a strong spin- 
up is measured within one 35 d cycle. Spectral characteristics observed during the X-ray dips are consistent with their interpretaion as 
due to partial covering as has been reported by several authors. The fundamental cyclotron absorption line is firmly observed in both 
pulse-averaged and pulse-phase resolved X-ray spectra. The energy, width, and the depth of the line are found to vary significantly 
with pulse phase. 

Key words. X-ray binaries ~ accretion disks - neutron stars 



1. Introduction 



Discovered in 1972 by the Uhuru satellite jGiacconi et al 
lTananbaumetalJ ll972') Her X-1 is one of the most intensively 
studied accreting pulsars. Being part of a low mass X-ray bi- 
nary it shows strong variability on very different time scales: 
the 1.24 s spin period of the neutron star, the 1.7 d binary pe- 
riod, the 35 d period of precessiori of the warped and tilted 
accretion disk jGerend & Bovntonlll976l: IShakuraet al] 119991: 
iKetsaris et al J2000h. and the 1 .65 d period of th e pre-ecUpse dips 
jGiacconi etalJll973l: iTananbaum et al.lll972h . Due to the high 
orbital inclination of the system (/ > 80°) the counter-orbitally 
precessing warped accretion disk around the neutron star cov- 
ers the X-ray source from the observer during a substantial part 
of the 35 d period. This gives rise to the alternation of so-called 
on (high X-ray flux) and ojf (low X-ray flux) states. The 35 d 
cycle contains two on states - the main-on and the short-on - 
separated by ~ 7 - 8 d ojf state. The X-ray flux in the mid- 
dle of the main-on state is ~ 4-5 times higher than that in 
the middle of the short-on. The sharp transition from the off- 
state to the main-on is called the turn-on of the source. Turn- 
ons are usually used for counting the cycles. The 35 d period 
manifest s itself also in variations of the shape of X-ray pulse 
profiles (ISoongetal.lll990bl: iTrumper et all 119861: iDeeteretalJ 
[1998"; 'Scott et al."2000') and modulat ion of opticaU ight curves 
(Gerend & Bovnton 1976; Howarth & Wilsonll 19831) . 



we present X-ray observations of Her X-1 performed by the 
INTEGRAL observatory. 

In Sect. |2] we describe the observations and data process- 
ing. In Sect. 13.11 we present energy- and time-dependent pulse 
profiles. In Sect. 13.21 we explore the behavior of the pulse pe- 
riod. The spectral analysis is presented in Sects. 14.11 (pulse- 
averaged spectra), 14.21 (spectral changes during X-ray dips), 
and 14.31 (pulse-resolved analysis). The results are discussed in 
Sect.|5] In Sect. |6] we present a summary and conclusions. 



2. Observations and data processing 



In spite of a large amount of available observational data, 
the physical interpretation of several observed phenomena in 
Her X-1 are still unclear. For example, the physical mechanisms 
responsible for the disk precession, the X-ray dips, and the evo- 
lution of X-ray pulse profiles are highly debated. In this work 



Her X-1 was observed by INTEGRAL dWinkler et al.l |2003') on 
July 22 - August 3, 2005 (MJD: 53573-53585 Q)- The obser- 
vations were spread over ~5 orbital revolutions of the satel- 
lite: revs. 338, 339, 340, 341, and 342. In ou r analysis we use d 
data obtained with the instruments JEM-X (Lund et al.ll2003h . 
IBIS/ISGRI (Lebrun etalJ|2003|; lUbertini et al.1 12003 1). and SPI 
( Vedrenne et al. 2003). We used the energy range ~3-20 keV 
for JEM-X and -20-100 keV for ISGRI and SPI (the spectrum 
of Her X-1 is falling off steeply with energy so that there is no 
signal above 100 keV). The data processing was performed with 
the version 6.0 Offline Science Analysis (OSA) software dis- 
tributed by ISDC (Courvoisier et al. 2003). An additional gain 
coiTection based on the analysis of the position of the tungste n 
background line was added (similar to iTsygankov et al]|2006l) . 
For the pulse-phase resolved spectroscopy we also used the soft- 
ware developed at lASF Palermo for IBISASGRE dMineo et al.1 



MJD = JD - 2400000.5 

,http : //www . ifc . inaf . it/~f errigno/INTEGRALsoftware .html | 
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Fig. 1. X-ray light curves of Her X-1 obtained with ISGRI (top 
panel) and JEM-X (bottom panel). The solid curve in the bottom 
panel show s the ASM RXTE lig ht curve averaged over many 
35 d cycles jKlochkov et al.ll2006[) . The INTEGRAL revolution 
numbers are indicated. 



120061: iFerrigno et"ani2007h . This software allows to construct 
phase-energy matrices from which pulse-resolved spectra and 
energy-resolved pulse profiles can easily be extracted. The re- 
sults of the analysis of pulse phase averaged spectra were in- 
dependently checked using the software package developed at 
Space Research Institute, Moscow. 

Figure [1] shows the light curves of Her X-1 obtained with 
ISGRI (top) and JEM-X (bottom). The observations partially 



cover a main-on state of the source, 35 d phases 



-0-0.11 



and 0pie ~0. 20-0. 28 (0pie = is the phase of the turn-on of the 
source). The turn-on occurred at MJD ~ 53574.7, corresponding 
to orbital phase 0orb ~ 0.7. The s ohd curve on top o f the JEM-X 
light curve is the ASM RXTE (iLevine et all 1 1 9961) light curve 
averaged over many 35 d cycles and renormolized to match the 
one from JEM-X. All JEM-X data from revolution 341 are re- 
jected by the ISDC team because of a strong solar flare. It is 
seen that the light curve is quite typical for Her X-1. Pre-eclipse 
dips (decrease of the flux just before the eclipse) are clearly seen. 



3. Timing analysis 

3.1. Pulse profiles 

To construct the 1.24 s pulse profiles of the source all photon 
arrival times were converted to the solar system barycenter and 
corrected for binary motion. The orbital parameters (Table [U 
used for the binary correction are based on an updated ephemeris 
arrived at by combining historical timing data with the timing 
results of the latest RXTE observations of Her X-1 (the cor- 
responding paper is being prepared). As a folding period for 
constructing pulse profiles we used Pspi„ - 1.23775836(10) s 
which was obtain ed with the standard epoch folding method 
(iLeahv et al.lll983h . 

It was found that the pulse profiles do not change signifi- 
cantly during the start (revs. 339 and 340) and during the end 
(rev. 341) of the main-on state. Therefore, we constructed the 



Table 1. Orbital parameters of Her X- 1 used to correct the arrival 
times of photons for orbital motion in the binary. 

Poyb = 1.700167233 d, 

= 53571.982111 (MJD), 
asin/ = 13.1831 lightsec. 



pulse profiles separately for these two intervals without fur- 
ther splitting in time. We excluded the data falling inside X- 
ray eclipses, pre-eclipse dips or located outside the main-on 
state (during the off'-state the data are very noisy and show only 
marginal sine-like pulsations). The resulting energy-resolved 
pulse profiles are shown in Fig. |2l It is seen that the profiles 
change significantly from the start to the end of the main- 
on state. The main peak (pulse phase ~0.6-0.8) is stronger 
and narrower at the end of the main-on. The interpulse (pulse 
phase -0.1-0.2) almost disappears. This behavior is typical for 
Her X-1: analogous pu lse profiles were o bserved by RXTE at 
different times (see e.g. lKuster et al.ll2005l) . We also studied the 
energy dependence of the pulse fraction F determined as 



(1) 



where /^ax and I^in are background-corrected count rates in the 
maximum and minimum of the pulse profile respectively. It is 
clearly seen (Fig. [3]l that the pulse fraction increases with en- 
ergy (there may be a saturation at energies above ~30 keV, how- 
ever the statistics is marginal). The structure of the profile is also 
energy-dependent, going from a triple-peaked shape at lower en- 
ergies to an almost single-peaked at higher energies. The main 
peak becomes narrower toward higher energies. A deeper anal- 
ysis of the pulse profiles in cluding the modeling of thei r shape 
(similar to that performed in Panchenko & Postnovll994l but tak- 
ing into account 35 d variation) is ongoing and will be presented 
elsewhere. 



3.2. Pulse period behavior 

To explore more precisely the intrinsic (not affected by the 
orbital motion in the binary) spin period of the neutron 
star Pspin during the INTEGRAL observations we performed 
a phase-con n ectio n analysis similar to Ferrigno et al. (200"^; 
iDeeter et al.l (OMlD using well defined average pulse profiles. 
The length of the time intervals used to produce a single pulse 
profile was 3-5 INTEGRAL Science Windows (6-10 ksec). This 
allowed to produce pulse profiles of a quality allowing to mea- 
sure their shift in pulse phase with respect to other pulse profiles 
within ~0.01 s uncertainties. Figure |4] shows the observed pulse 
arrival times minus the calculated times (using a constant pulse 
period). If we assign the same statistical uncertainty of 0.01s (de- 
termined from the scattering of data points) to each data point 
then the linear fit to the data corresponding to a constant pulse 
period (dashed line) gives x^^^ - 1-9 for 23 d.o.f. while the 
quadratic fit corresponding to the presence of a non-zero Pspin 
results inx^^^ - l.Q for 22 d.o.f. We conclude, therefore, that the 
neutron star is spinning up with -/"spin = (5.8 + 1.5) x lO^'-'s/s. 
This value is significantly higher than the averaged spin-up trend 
of Her X-1 which is ~ 1.1 x 10"'^ s/s. Around the average spin- 
up strong variability of the pulse period from one main-on to the 
next is known to exist (see e.g. Fig. 1 in lStaubert et al.ll2006l) . 
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Fig. 2. Energy-resolved X-ray pulse profiles of Her X-1 at the start of the main-on (left column) and at the end of the main-on (right 
column). 



4. Spectral analysis 

The spectral analysis has been perforrned u sing the XSPEC 
V. 1 1 .3.21 spectral fitting package ( lArnaudl ll996). For the start of 
the main-on the data from all three X-ray instruments have been 
used. The fit, however, is mainly driven by JEM-X (< 20 keV) 
and ISGRI (> 20 keV) data. The specti-a obtained with SPI 
(> 20 keV) have much lower statistics. For the end of the main- 
on only the higher energy part of the spectrum is available since 
the JEM-X data could not be used due to a solar flare (see above). 
Following the OSA User Manualfl we added systematic errors 
at a level of 2% in quadrature to all JEM-X and ISGIU spec- 
tral points in order to account for uncertainties in the response 
matrices of the respective instruments. 

4.1. Pulse-averaged spectra 

For constructing the pulse-averaged spectrum, we used the data 
from revolutions 339 and 340 (the start of the main-on state, see 
Fig.[T]i. During this interval the data from all three X-ray instru- 
ments are available. Eclipses and dips, which are clearly seen 



|http : //isdc . unige . ch/index . cgi?Support+documents| 



in the light curve (Fig. [U and have sharp boundaries, were ex- 
cluded from the analysis. To model the broad band spectral con- 
tinuum we used a power law with the exponential cutoff model 
highecut of XSPEC: 



E^^ ■ exp ^- 



if £■ < /icutoff 

if £■ > -Ecutoif, 



(2) 



where E is the photon energy; F, iScutoff, and isfoid are model pa- 
rameters. A Gaussian line is added to model the iron fluores- 
cence line at ~6.5 keV. The cyclotron absorption feature is mod- 
eled by a multiplicative absorption line with a Gaussian optical 
depth profile. So, the final spectral function / is the following: 



■ exp < -Tcyci exp 



(E £'cycl) 
^^cycl 



2^ 



+ A" exp 



(E - £Fe ) 



2^ 



(3) 



where E^yd and cTcyd are the centroid energy and width of the 
cyclotron line, fipe and crpe are the energy and width of the iron 
emission line. T^yci and K are the numerical constants describing 
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Fig. 3. Pulse fraction of Her X-1 as a function of energy. 



Table 2. Best fit spectral parameters of Her X-1 for the observa- 
tions of revs. 339 and 340. lcr(68%)-uncertainties {x^i„ + 1) for 
one parameter of interest are shown. 



Parameter 



^cutoff [keV] 
£:foid [keV] 
^Fe [keV] 
a-Fe [keV] 

^cycl 

F JEM-X 
F ISGRI 
^SPI 



Value 



0.91 ±0.01 

25 5+"-2 

'•"-0.2 

6 57+" " 

"•-"-0.13 

0.50!°- 
38.2!«;« 

Q ^+0.4 

1.0 (fixed) 
0.8lt»^ 
0.96 ± 0.02 
1.19/207 



0.04 



0.02- 




53576 53578 53580 53582 
Time [MJD] 

Fig. 4. Observed pulse aiTival times of Her X- 1 minus calculated 
times using a constant pulse period. The dashed line shows a 
linear fit corresponding to a constant pulse period while the solid 
line shows a quadratic fit corresponding to a /'spin - (-5.8+1.5)x 

lO-'^s/s . 



the strength of the cyclotron line and iron line, respectively. The 
count rate spectrum is shown in Fig.|5] The middle panel shows 
the residuals after fitting the spectrum with the continuum model 
without the iron emission line and the cyclotron absorption line. 
Systematic features around 6 keV and 40 keV are clearly seen. 
The bottom panel shows the residuals of our final fit which in- 
cludes the two lines. 

To account for large systematic uncertainties in the nor- 
malization of the instruments we introduced in our models a 
free multiplicative factor for each instrument: Fisgri, ^^spi and 
^^jEM-x (for JEM-X the factor was fixed to 1.0). The best-fit pa- 
rameters with corresponding lcr(68%)-uncertainties are listed in 
Table |2] The position of the cyclotron line at 38.2+Qg keV is 
consistent with that measured by RXTE at a similar time (see 
Table 1 in Staubert et al. 2007b) and, therefore, supports the cor- 
relation between the cycl otron line energy and the maximum 
main-on flux found by Staubert et al.l (i2007b) . 

We have checked the presence and the centroid energy of the 
cyclotron line using different continuum models (such as e.g. 
Fermi-Dirac cutoff, Tanaka 1986) as well as using different line 
profiles (e.g. Lorenzian profile). It was found that the presence 
and the energy of the feature are independent of the choice of 



the continuum and the line profile model. The highecut model 
was found to provide a better description of the spectral con- 
tinuum in Her X-1 with respect to other models. Additionally, 
it allows to compare our results with previous observations of 
the source, mo st of which were a nalyzed using the highecut 
model (see e.g. iGruber et al.ll2001i) . A discontinuity of the first 
derivative of the spectral function at ^cutoff ~ 25 keV (which can 
be noticed in the residuals in Fig. |5]l is far from the cyclotron 
line energy and does no t affect the line parameters (see however 
iKretschrnar et al.lll997b . 

We also explored the pulse-averaged spectrum obtained dur- 
ing revolution 341 (end of the main-on). The continuum spectral 
parameters are poorly constrained in this case due to the absence 
of the low-energy (< 20keV) part of the spectrum. However, the 
found position of the cyclotron line Ecyci - 37.3^} j keV is con- 
sistent with that measured in revolutions 339 and 340. 

4.2. X-ray dips 

In the X-ray light curve taken with INTEGRAL one can see 
pre-eclipse dips on three subsequent orbits after the turn-on 
(Fig. IB- The turn-on of the source can be considered as egress 
fr om an anomalous dip at (po^b ~ 0.5 (according to the model 
of' Shakuraet"aLl (1 19991) the source is partially screened by the 
outer accretion disk rim in both cases: during the turn-on of the 
source and during anomalous dips). Thus, the light curve ob- 
tained with INTEGRAL contains three pre-eclipse dips and one 
anomalous dip. Using the JEM-X data we calculated the ratio 
H = /8-20kev//3-8keV, whcrc /8-2()keV and /3-8keV are the count 
rates in the respective energy ranges. The ratio // as a function 
of time is shown in the bottom panel of Fig. |6] It can be seen 
that H increases during X-ray dips most probably indicating the 
low-energy absorption by the cold matter in the accretion stream 
(during the pre-eclipse dips) and in the outer rim of the disk (dur- 
ing the anomalous dip). Even though the result is dominated by 
the last two pre-eclipse dips, the data of the first pre-eclipse dip 
and the anomalous dip are consistent. 

To check the absorption hypothesis we constructed JEM-X 
and ISGRI spectra corresponding to the four dips and compared 
them with those obtained outside the dips. We are aware of the 
fact that due to different formation mechanisms of anomalous 
and pre-eclipse dips their spectral characteristics might be dif- 
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Fig. 5. Top: The pulse-averaged X-ray spectrum of Her X-1 ob- 
tained during revolutions 339 and 340 (start of the main-on); 
middle: residuals after fitting the spectrum with the continuum 
model described by Eq. |2] without the iron emission line and 
cyclotron absorption line (to avoid confusion the SPI data are 
not shown in this panel, their statistics is lower than that of 
IBIS/ISGRI and JEM-X); bottom: the residuals after including 
the iron and cyclotron lines in the model. 
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Fig. 6. The JEM-X 3-20 keV light curve ( top panel) and the ratio 
of count rates in the harder (8-20 keV) and softer (3-8 keV) 
energy ranges as a function of time (bottom panel). 



ferent. But due to the low statistics we cannot perform a separate 
analysis of the two kinds of dips (the spectral parameters are not 
restricted in this case). To fit the spectrum from the dips we used 
the approach of Kuster et al. (2005). It was assumed that a com- 
bination of direct and absorbed radiation is observed during the 
dips. Thus, we used a partial covering model which combines 



both, absorbed and non-absorbed spectra to fit the data during 
the dips. We used the same continuum model as in the previous 
section. The final spectral model can be written as 

/(£)=A-[l+a(£)-W], (4) 

where /com is the continuum model described in the previous 
section, and 

a(E) = B ■ e-^H^r„(£)^ 

where crbf(£') is the photoabsorption cross-section 
per hydrogen atom for matter of cosmic abundances 
(iBalucihska-Church & McCammonI [T992I) used in the phabs 
model of XSPEC and A^h is the equivalent hydrogen column 
density. A larger value of B implies a larger degree of ab- 
sorbed flux. From the fit of the spectrum of the dips we found 
B = 0.24 + 0.03 and A^h = (111![9) x 10^2 cm^^. The best-fit 
parameters of the spectral continuum /cont are consistent with 
those found from fitting the spectrum outside the dips (Table|2]i. 
To compare this result with the spectrum outside the dips, the 
latter was also fit by the partial covering model (Eq. |4]i. Both 
spectra (inside and outside the dips) are shown in Fig. [7] For the 
spectrum outside the dips (the upper one) the A^h parameter was 
found to be consistent with zero and the B factor, consequently, 
was not restricted. We conclude therefore, that the X-ray data 
of Her X-1 during the dips are consistent with the partial 
absorption model. 

4.3. Pulse-phase-resolved spectra 

It is well known th at the X-ray spectrum of Her X-1 varies with 
1.24 s pulse phase dVoges et alJI 19821: IVogeslll985HSoong et all 
^19903^. So, we have performed a separate analysis of the spectra 
accumulated in different pulse phase intervals. Figure |8] shows 
an example of pulse resolved spectra of the source. Variability 
of the continuum and the cyclotron line is clearly seen. 

Since the shape of the pulse profile is changing significantly 
from the start of the main-on (revs. 339, 340) to its end (rev. 341), 
we analyzed the data from these two intervals separately. Phase 
binning in each case was chosen to provide a similar statistics in 
each spectral bin. Pulse phase zero is the same as in Fig.|2] The 
X-ray spectrum of each phase bin was fitted with the spectral 
model described in Sect. 14.11 Figures |9] and [TOl show the best-fit 
spectral parameters as a function of pulse phase for revs. 339 and 

340 (left) and the rev. 341 (right). Vertical eiTorbars coiTespond 
to lcr(68%)-uncertainties. The dotted line shows the coiTespond- 
ing pulse profile. The spectral parameters in different pulse phase 
bins are also listed in Table[3] Since no JEM-X data are available 
for rev. 341 the power law photon index F could not be restricted 
by the fit in this case. The iron line at ~6.5 keV was detected only 
in the phase bin 0.1-0.54, corresponding to the off-pulse inter- 
val. The data from other phase intervals do not require inclusion 
of the line in the spectral model to obtain a good fit. For the data 
coiTesponding to revs. 339 and 340, it is seen that all continuum 
and the cyclotron line parameters are highly variable with pulse 
phase. For revolution 341 the picture is not so clear due to the 
absence of the low-energy (< 20 keV) part of the spectrum. Also 
the total exposure time and the average flux of the source in rev. 

341 are smaller than those during revolutions 339 and 340. 

5. Discussion 

5. 1 . Pulse profiles and pulse period 

In Sect. 13.11 we constructed 1.24 s X-ray pulse profiles of 
Her X-1. The shape of the profiles is both energy- and time- 
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Table 3. Best fit spectral parameters of Her X- 1 in different pulse phase intervals. 1 cr(68%)-uncertainties 
of interest are shown. 
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Fig. 7. Examples of pulse-averaged spectra of Her X-1 from the 
dips (lower curve) and outside the dips (upper curve). The dif- 
ference of the spectral shape at lower energies caused by photo- 
absorption inside the dips is clearly seen. The step-like feature 
at ~7 keV in the spectrum of the dips corresponds to the iron 
K-edge (as implemented in the phabs model of XSPEC). 



dependent. At higher energies the main peak becomes narrower 
and the pulsed fraction increases (Figs.|2]and[3ll. Such a depen- 
dence o f the profile on energy is typical for accreting pulsars 
(see e.g. iTsveankov et aL I l2007h . It can be understood in a sim- 
ple purely geometrical picture: if the rotation axis of the neutron 
star is inclined with respect to the axis of its magnetic field, it 
is expected that the upper part of the accretion column emitting 
softer photons is seen during a larger part of the neutron star spin 
period while the emission region of harder photons, the "base" 
of the accretion column, is screened by the neutron star surface 
during most of the spin period. Additionally, the emission dia- 
gram of harder photons most probably originating close to the 
base of the accretion column is believed to be narrower than that 
of softer photons due to the strong dependence of the scattering 
cross-section on the angle between the direction of a photon and 
the magnetic field lines (see e.g. Basko & Sunvaev 1976). 

The difference of the pulse profiles in the left and right 
columns of Fig. |2] (corresponding to the start and the end of 




10 20 
channel energy (keV) 



Fig. 8. Examples of pulse-resolved X-ray spectra of Her X-1 
during revolutions 339 and 340. Pulse phase intervals are indi- 
cated (the pulse phase zero is the same as in Fig.|2|i. 



the main-on state, respectively) shows the evolution of the pro- 
files with 35 d pha se. Two models have been proposed to explain 
these variations. Scott et al.l (1200 0) suggested that the change 
of the profile from the start to the end of the main-on state is 
due to a "resolved occultation" of the emitting regions on the 
neutron star by the precessing accretion disk. According to this 
model, the disk progressively occults the neutron star towards 
the end of the main-on. The other mo del assumes free preces- 
sion ofJhejieutron_st^J Trum 1986; Kahabka 198^ 
ll989l:lProkhorov et al.lll990tlShakura et al.lll998l) . In this model 
the observed behavior of the profile is explained assuming that 
the emission region on the star surface has a complex shape 
(due to the presence of higher multipole components of the neu- 
tron star's magnetic field). Changing the viewing conditions of 
the emitting region with the phase of the fr ee precession causes 
the observed varia tion of the pulse profile (iKetsaris et al.ll2000l; 
IWilms et alj|2003l) . In pulse profiles obtained with INTEGRAL 
one can see that the count rate at the maximum of the main peak 
is higher at the end of the main-on than that during the start of 
the main-on (Fig.|2]). In the case of a progressive occultation one 
would only expect a decrease in the intensity of any feature in 
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the profile towards the end of the main-on. In the model of free 
precession, however, viewing conditions of different parts of the 
emission region may change with the precessional phase of the 
star in such a way that the intensity of particular features of the 
profile (e.g of the main peak) will increase towards the end of the 
main-on. We argue, therefore, that the INTEGRAL observations 
analyzed in this work support the model of a freely precessing 
neutron star as an explanation of the time variation of the pulse 
profile in Her X- 1 . 

In Sect. 13.21 we determined the time derivative of the in- 
trinsic (not affected by the orbital motion in the binary) pulse 
period of the neutron star in Her X-1 during the INTEGRAL 
observation. The pulsar was found to spin-up with the rate 
-^'spin = (5. 8 + 1. 5)x lO"''' s/s. This value is more than five times 
larger than the mean spin-up trend of Her X-1 (~ 1.1 x 10"'^' s/s) 
and consistent with the variations of Pjpin between subsequent 
35 d cycles, e g. as seen by BATSE (Kunz 1996; Nagaselll989t 
ISunvaev et al.ll988l) and RXTE ( Staubert et al.l2006 >). However, 
this is the first time that such strong pulse period variations have 
been measured within one 35 d cycle. A deeper analysis of these 
pulse period variations and their corelation with the X-ray lu- 
minosity (including the data from INTEGRAL presented here) 
is currently being performed and will be presented in a separate 
paper. 



5.2. Absorption during X-ray dips 

In Sect. 14.21 we have shown that the X-ray spectrum of Her X- 1 
obtained during X-ray dips can be modeled using a par- 
tial covering model which assumes that the observed spec- 
trum is a combination of direct and absorbed radiation. Such 
a combination has already been observed in Her X-1 during 
X-ray turn-ons (Q Pavison & Fabian 1977; iBecker et al.lll977l: 
iParmar e t al. 1980; Kuster et al. 120051) and also during X-ray 
dips IVrtilek & Halpern ,1985. although the observations in- 
cluded only the lower energy part of the spectrum <20 keV). 

According to the model of IShakura et al.l ( 11999 ') and 
iKlochkov et al.l (l2006l) (which we adopt in this work), most of 
the X-ray dips are produced by the occultation of the X-ray 
source by the cold matter in the accretion stream which moves 
out of the system's orbital plane. If the stream is not homoge- 
neous but consists of "blobs" of material then the source will be 
screened during some time intervals within a dip. Between these 
intervals the source will be visible or seen through less dense 
material. This picture seems to be confirmed by the complicated 
behavior of the light curve inside the dips (Fig.|6]l. In our analy- 
sis, in order to construct an X-ray spectrum of acceptable quality 
we had to accumulate the flux during complete dips. The spec- 
trum in this case will contain a superposition of the absorbed and 
non-absorbed flux. The partial covering model fits the spectrum 
of the dips clearly better than the simple absorption model. 



5.3. Pulse-pliase variabiiity of tine X-ray spectrum 

The dependence of the spectral parameters in Her X- 1 on pulse 
phase observed by INTEGRAL bas ically confirm s the results 
obtained earUer with HEAO-1 ( Soong et alJll990bl) but provide 
better restriction of continuum parameters. The spectral varia- 
tion with pulse phase is a common feature in X-ray pulsars (see 
e.g. Kreykenbohm et al. 2004, and references therein) which is 
usually attributed to the change of the viewing angle of the ac- 
cretion region on the neutron star surface. Below we show that 
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Fig. 9. Best-fit spectral continuum parameters as functions of the 
pulse phase for revolutions 339 and 340 (left) and the revolution 
341 (right). The dashed curve shows the coiTesponding pulse 
profile. Vertical eiTor bars coiTespondto lcr(68%)-uncertainties. 



the pulse-phase variability observed in Her X-1 is basically con- 
sistent with this interpretation. 

As one can see in Fig. |9] the power law photon index F de- 
creases during the main peak. It ranges from ~1.0 in the off- 
pulse to ~0.5 close to the maximum of the main peak (see also 
ILutovinov et al.ll2000l) . This eff'ect reflects the sharpening of the 
main peak with energy which is observed in energy-resolved 
pulse profiles constructed in Sect. l3.1l (Fig.l2]i. The spectral hard- 
ening in the main peak can be explained by the dependence of 
the optical depth on t he angle between th e line-of-sight and the 
magnetic field lines dPravdo et al.lll97'7l) . The closer the view- 
ing direction is to the magnetic axis, the deeper we look into 
the emission region where harder photons originate. Since it is 
generally accepted that the main peak coiTesponds to the radia- 
tion from one of the two magnetic poles, the viewing direction is 
closest to the magnetic axis during the main peak which causes 
the observed effect. 

The exponential folding energy Em^ is much higher during 
the main peak than during the rest of the pulsation cycle. In 
comptonized X-ray spectra (which are normally observed from 
accreting pulsars) this parameter is believed to be proportiona l 
to the plasma temperature (see e.g.'Rvbic ki & Lightmanll 19791) . 
As already mentioned, during the main peak we see the emission 
from higher optical depth. This means that the observed radia- 
tion originates closer to the base of the accretion column where 
the plasma temperature is higher. 

The cyclotron line centroid energy E^yd is also higher dur- 
ing the main peak of the profile (Fig. [TO] i. A relative amplitude 
of its variation is ~25%. As argued from F and Em^, during the 
main peak we see X-rays that originate closer to the neutron star 
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Fig. 10. The same as in Fig.|9]but for the cyclotron hne parame- 
ters (left: revolutions 339 and 340; right: revolution 341). 



surface, i.e. where the magnetic field strength is higher. The in- 
crease in Ecyc\ during the main peak is, therefore, qualitatively 
consistent with this picture. If one assumes that the radiation 
comes from a compact accretion column located at the magnetic 
pole of the neutron star with a pure dipole field then the change 
in the height of the observed emission region necessary to pro- 
duce a ~25% variation of the magnetic field strength is km 
(for a neutron s tar radi us of 10 km). However, as it was shown by 
IStaubert et al.l (l2007al) . the height of the emission and line form- 
ing region in case of Her X-1 is most probably much smaller, 
~ 10"* cm. Thus, a changing height of the emitting region above 
the neutron star cannot explain the observed variability of the 
cyclotron line energy with pulse phase in this source. 

Another possibility for the line energy variability is to as- 
sume a complicated structure of the magnetic fiel d at th e site 
of X-ray emission. As show n by Shakura et al.l (1199 Ih and 
iPanchenko & Postnovl d 19941) . the complex shape of the ob- 
served pulse profiles in Her X- 1 suggests that such a complicated 
field structure including higher multipole components is indeed 
present in the source. In this case, at different rotational phases 
of the neutron star we will observe emission coming from the re- 
gions corresponding to different sub-structures of the non-dipole 
magnetic field with very different field strengths. Furthermore, 
the numerical modeling of Her X-I pulse pr ofiles observed with 
EXOSAT performed by IPanchenko & Postnovl (1 19941) suggests 
that the magnetic field stength at the ring-like structure on the 
neutron star's surface corresponding to the quadrupole compo- 
nent of the field (which was introduced in the model) is indeed 
by ~27% lower than that at the magnetic pole. Thus, a compli- 
cated magnetic field structure on the surface of the neutron star 
can cause the observed 25%-variation of E^yd with pulse phase. 



6. Summary and conclusion 

In this work we analyzed the INTEGRAL observations of the ac- 
creting X-ray pulsar Her X-1. X-ray pulse profiles of the source 
are constructed and their energy- and time-dependence are dis- 
cussed. A strong spin-up is found during the INTEGRAL obser- 
vations. This is the first time that pulse period variation is mea- 
sured within one main-on state of the source. The value of the 
spin-up rate is consistent with the typical difference between the 
values of the pulse period previously found in adjacent 35 d cy- 
cles. Spectral changes during X-ray dips are studied. The X-ray 
spectrum of the dips was modeled by a partial covering model 
which assumes that the observed spectrum is a combination of 
direct and absorbed radiation. This is in agreement with the 
previously suggested interpretation assuming that the obscuring 
matter (accretion stream) is n ot homogeneous but rathe r consists 
of blobs of material (see e.g. lVrtilek & HalpernllI985l) . Energy- 
resolved X-ray pulse profiles as well as the variation of cyclotron 
line and continuum parameters with 1 .24 s pulse phase were ex- 
plored. The spectral changes with pulse phase are shown to be 
qualitatively (for the cyclotron line energy - also quantitatively) 
consistent with those expected from the viewing conditions of 
the complex emitting region which vary with the rotational phase 
of the neutron star. 

Acknowledgements. This research is based on observations with INTEGRAL, 
an ESA project with instruments and science data centre funded by ESA 
member states (especially the PI countries: Denmark, France, Germany, Italy, 
Switzerland, Spain), Czech Republic and Poland, and with the participation of 
Russia and the USA. The work was supported by the DFG grants Sta 173/31-2 
and 436 RUS 1 13/717/0-1 and the coiTesponding RBFR grants RFFI-NNIO-03- 
02-04003 and RFFI 06-02-16025, as well as DLR grant 50 OR 0302 and RFFI 
07-02-01051. We also thank ISSI (Bern, Switzerland) for its hospitality during 
the team meetings of our collaboration. AL acknowledges the financial support 
of the Russian Science Support Foundation. 



References 

Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference Series, 

Vol. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. 

lacoby & J. Barnes, 17 
Balucihska-Church, M. & McCammon, D. 1992, ApJ, 400, 699 
Basko, M. M. & Sunyaev, R. A. 1976, MNRAS, 175, 395 
Becker, R. H., Boldt, E. A., Holt, S. S., et al. 1977, ApJ, 214, 879 
Courvoisier, T. J.-L., Walter, R., Beckmann, V., et al. 2003, A&A, 411, L53 
Davison, P J. N. & Fabian, A. C. 1977, MNRAS, 178, IP 
Deeter, J. E., Pravdo, S. H., & Boynton, P E. 1981, ApJ, 247, 1003 
Deeter, J. E., Scott, D. M., Boynton, P E., et al. 1998, ApJ, 502, 802 
Ferrigno, C, Segreto, A., Santangelo, A., et al. 2007, A&A, 462, 995 
Gerend, D. & Boynton, P E. 1976, ApJ, 209, 562 
Giacconi, R., Gursky, H., Kellogg, E., et al. 1973, ApJ, 184, 227 
Gruber, D. E., Heindl, W. A., Rothschild, R. E., et al. 2001, ApJ, 562, 499 
Howarth, I. D. & Wilson, B. 1983, MNRAS, 202, 347 
Kahabka, P 1987, NASA STI/Recon Technical Report N, 88, 19405 
Kahabka, P 1989, in ESA SP-296: Two Topics in X-Ray Astronomy, Volume 1: 

X Ray Binaries. Volume 2: AGN and the X Ray Background, ed. J. Hunt & 

B. Battrick, 447^52 
Ketsaris, N. A., Kuster, M., Postnov, K., et al. 2000, in Proc. Int. Workshop "Hot 

Points in Astrophysics", JINR, Dubna, p. 192, ed. V. Belyaev, arXiv:astro- 

ph/0010035 

Klochkov, D. K., Shakura, N. I., Postnov, K. A., et al. 2006, Astronomy Letters, 
32, 804 

Kretschmar, P., Kreykenbohm, I., Wilms, J., et al. 1997, in ESA Special 
Publication, Vol. 382, The Transparent Universe, ed. C. Winkler, T. J.-L. 
Courvoisier, & P. Durouchoux, 141 
Kreykenbohm, I., Wilms, J., Coburn, W., et al. 2004, A&A, 427, 975 
Kunz, M. 1996, PhD thesis. University of Tubingen, Germany 
Kuster, M., Wilms, J., Staubert, R., et al. 2005, A&A, 443, 753 
Leahy, D. A., Eisner, R. F, & Weisskopf, M. C. 1983, ApJ, 272, 256 
Lebrun, F, Leray, J. P, Lavocat, P, et al. 2003, A&A, 41 1, L141 
Levine, A. M., Bradt, H., Cui, W., et al. 1996, ApJ, 469, L33 
Lund, N., Budtz-J0rgensen, C, Westergaard, N. J., et al. 2003, A&A, 411, L231 



D. Klochkov et al.: INTEGRAL observations of Her X-1 



9 



Lutovinov, A. A., Grebenev, S. A., Pavlinsky, M. N., & Sunyaev, R. A. 2000, 

Astronomy Letters, 26, 691 
Mineo, T., Ferrigno, C, Foschini, L., et al. 2006, A&A, 450, 617 

Nagase, F 1989, PASJ, 41, 1 

Panchenko, 1. E. & Postnov, K. A. 1994, A&A, 286, 497 
Parmar, A. N., Sanford, P W., & Fabian, A. C. 1980, MNRAS, 192, 311 
Pravdo, S. H., Boldt, E. A., Holt, S. S., & Serlemitsos, R J. 1977, ApJ, 216, L23 
Prokhorov, M. E., Shakura, N. I., & Postnov, K. A. 1990, Thirty-five-day cycle of 

HER X-1 : Synthesised optical light curves in the model of freely processing 

neutron star. Tech. rep. 
Rybicki, G. B. & Lightman, A. P. 1979, Radiative processes in astrophysics (New 

York, Wiley-lnterscience, 1979. 393 p.) 
Scott, D. M., Leahy, D. A., & Wilson, R. B. 2000, ApJ, 539, 392 
Shakura, N. 1., Postnov, K. A., & Prokhorov, M. E. 1991, Soviet Astronomy 

Letters, 17, 339 

Shakura, N. 1., Postnov, K. A., & Prokhorov, M. E. 1998, A&A, 331, L37 
Shakura, N. 1., Prokhorov, M. E., Postnov, K. A., & Ketsaris, N. A. 1999, A%A, 

348,917 

Soong, Y., Gruber, D. E., Peterson, L. E., & Rothschild, R. E. 1990a, ApJ, 348, 
641 

Soong, Y, Gruber, D. E., Peterson, L. E., & Rothschild, R. E. 1990b, ApJ, 348, 

634 

Staubert, R., Klochkov, D., & Rodina, L. 2007a, A&A, [in preparation] 
Staubert, R., Schandl, S., Klochkov, D., et al. 2006, in American Institute 
of Physics Conference Series, Vol. 840, The Transient Milky Way: A 
Perspective for MIRAX, ed. J. Braga, F. D' Amico, & R. E. Rothschild, 65-70 
Staubert, R., Shakura, N. 1., Postnov, K., et al. 2()07b, A&A, 465, L25 
Sunyaev, R. A., Gilfanov, M. R., Churazov, E. M., et al. 1988, Soviet Astronomy 
Letters, 14, 416 

Tanaka, Y. 1986, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 255, 
lAU Colloq. 89: Radiation Hydrodynamics in Stars and Compact Objects, ed. 
D. Mihalas & K.-H. A. Winkler, 198 
Tananbaum, H., Gursky, H., Kellogg, E. M., et al. 1972, ApJ, 174, L143 
Trumper, J., Kahabka, P, Oegelman, H., Pietsch, W., & Voges, W. 1986, ApJ, 
300, L63 

Tsygankov, S. S., Lutovinov, A. A., Churazov, E. M., & Sunyaev, R. A. 2006, 
MNRAS, 371, 19 

Tsygankov, S. S., Lutovinov, A. A., Churazov, E. M., & Sunyaev, R. A. 2007, 

Astronomy Letters, 33, 368 
Ubertini, P, Lebrun, F., Di Cocco, G., et al. 2003, A&A, 41 1, L131 
Vedrenne, G., Roques, J.-P, Schonfelder, V., et al. 2003, A&A, 411, L63 
Voges, W. 1985, NASA STl/Recon Technical Report N, 85, 34112 
Voges, W., Pietsch, W., Reppin, C, et al. 1982, ApJ, 263, 803 
Vrtilek, S. D. & Halpern, J. P 1985, ApJ, 296, 606 

Wilms, J., Ketsaris, N. A., Postnov, K. A., et al. 2003, Izvestiya Akademii Nauk, 

Ser. Fizicheskaya, 67, 310 
Winkler, C, Courvoisier, T. J.-L., Di Cocco, G., et al. 2003, A&A, 41 1, LI 



